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ABSTRACT 

The magnetospheres of accreting neutron stars develop electrostatic gaps with huge potential drops. Protons and 
ions, accelerated in these gaps along the dipolar magnetic field lines to energies greater than 100 TeV can impact 
onto the surrounding accretion disc. A proton-induced cascade develops, and charged pion decays produce v- 
emission. With extensive disc shower simulations using DPMJET and GEANT4, we have calculated the resulting 
//-spectrum. We show that the spectrum produced out of the proton beam is a power law. We use this result to 
propose accretion-powered X-ray binaries (with highly magnetized neutron stars) as a new population of point-like 
//-sources for km-scale detectors such as ICECUBE. As a particular example we discuss the case of A0535+26. 

We show that ICECUBE should find A0535+26 to be a periodic //-source, one for which the formation and loss 
of its accretion disc can be fully detected. Finally, we briefly comment on the possibility that smaller telescopes, 
like AMANDA, could also detect A0535+26 by folding observations with the orbital period. 

Subject headings: X-ray binaries: general, gamma-rays: theory, gamma-rays: observations, neutrinos: 
observations. X-ray binaries: A0535+26 

1. INTRODUCTION Compton Gamma Ray Observatory detected a 33-day, broad, 


X-ray binaries have fascinated those looking for Galactic 
neutrino sources (e.g. Berezinsky et al. 1985, Kolb et al. 1985, 
Gaisser and Stanev 1985). The basic idea is to use somehow 
the secondary object in the system to accelerate protons, which 
then could collide within a higher density medium. One pos¬ 
sibility that was earlier explored was to use the primary itself 
as the target for those accelerated hadrons (Berezinsky et al. 
1990). Depending on the grazing angle of the colliding pro¬ 
tons, on the size and type of the primary star, and on the effec¬ 
tiveness of the acceleration mechanism, 7 -rays and neutrinos 
(or us) could escape from the system. 7 -rays, however, could 
more naturally be produced in the accretion disc surrounding 
the neutron stars; detailed models for this possibility were pre¬ 
sented by Cheng and Ruderman (1989). Herein, we show that 
accreting X-ray binaries in which the compact object is a mag¬ 
netized neutron star are sources of high energy neutrinos that 
can be detected by forthcoming neutrino telescopes. Moreover, 
we show that the signal-to-noise ratio can be high enough as to 
allow timing studies and multiwavelength comparison. 

2. THE ACCRETION DISC OF A0535+26 

Of all the X-ray binaries, A0535+26 is one of the most stud¬ 
ied. A0535+26 is a Be/X-ray transient where the compact ob¬ 
ject is a 104 second pulsar in an eccentric orbit around the 
B0III star HDE 245770. Be stars are rapidly rotating objects 
that eject mass, irregularly forming gaseous discs in their equa¬ 
torial planes. Strong and recurrent X-ray outbursts were ob¬ 
served separated by 111 days, which has been identified with 
the orbital period (Giovannelli and Sabau Graziati 1992). These 
outbursts occur when the accretion onto the neutron star in¬ 
creases at periastron passage. (The average ratio of the X-ray 
luminosity at the periastron to that of apoastron is ~ 100 , 
Janot-Pacheco et al. 1987.) The BATSE instrument of the 


quasi-periodic oscillation in the power spectra of the X-ray flux, 
definitively showing that an accretion disc is formed during gi¬ 
ant outbursts (Finger et al. 1996). The 7 -ray light curve maxi¬ 
mum for this object is anti-correlated with X-rays at periastron 
passage (Hartman et al. 1999, Romero et al. 2001). To ex¬ 
plain this anti-correlation, a variation in the disc grammage is 
invoked: when the disc is fully formed, and the X-ray luminos¬ 
ity is at its maximum, the disc grammage is too high as to allow 
7 -ray photons above 100 MeV to escape. We now describe how 
the accretion disc model works. 

Accretion discs can penetrate the stellar magnetospheres of 
accreting rotating neutron stars (Ghosh and Lamb 1979). This 
penetration creates a broad transition zone between the unper¬ 
turbed disc flow —far from the star— and the co-rotating mag- 
netospheric flow —close to the star. In the transition zone, with 
inner radius ro, the angular velocity is Keplerian. Between r () 
and the co-rotation radius r co there is a thin boundary layer 
where the angular velocity significantly departs from the Keple¬ 
rian value. At r co the disc is disrupted by the magnetic pressure 
and the accreting mass is channelled by the field to impact onto 
the stellar surface, producing strong X-ray emission. 

The magnetosphere of an accreting neutron star spinning 
slower than that of the accretion disc can be divided into three 
regions: i) a region coupled to the star by the magnetic field 
lines that do not penetrate the disc, ii) an equatorial region 
linked to the disc by the field attached to it, and iii) a gap en¬ 
tirely empty of plasma separating both regions (Cheng and Ru¬ 
derman 1989). The magnetospheric region penetrated by field 
lines which do not intersect the disc co-rotates with the star at 
an angular velocity f1*. When the angular velocity of the disc 
0 ( i exceeds the angular velocity of the star (O* < 0 ( [), the equa¬ 
torial plasma between the inner accretion disc radius ry and the 
Alfven radius r,\ co-rotates with the disc to which it is linked 
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Fig. 1.— Sketch (not to scale) of the magnetosphere of an accreting magnetized neutron star, with the electrostatic gap produced when Q* < Protons entering 
into the gap from the region that co-rotates with the star (-) are accelerated along the field lines and impact onto the disc, initiating a shower, us and 7 -ray photons, 
resulting from pion decay, emerge from the opposite side of the disc and will generically be beamed due to momentum conservation. This figure is adapted from 
several of the works by Cheng and Ruderman, quoted in the text. 


by the frozen field lines. Inertial effects result in a charge sep¬ 
aration around the “null surface", 12* • B = 0. This leads to the 
formation of an electrostatic gap, with no charge (see Fig. 1). In 
this gap E ■ B t^O, and a strong potential drop is established. The 
maximum potential drop along the magnetic field lines through 
the gap is (Cheng and Ruderman 1989) 

Vmax ~ 4 X 10 14 /r 5/2 ' R^L^B ^ 1 V, (1) 

where M* ~ 1.4 —2.7 M 0 is the neutron star mass (Giovan- 
nelli and Sabau Graziati 1992), R(, = /?*/10 6 cm is its radius, 
B \2 = B, /10 12 G is the magnetic field of the star, L 37 is the X- 
ray luminosity in units of 10 37 erg s -1 , and (3 = 2ro/r A is twice 
the ratio between the inner accretion disc radius and the Alfven 
radius. The Alfven radius for spherical accretion can be deter¬ 
mined from the condition that the unscreened magnetic energy 
density of the stellar field becomes comparable to the kinetic 
energy density of the accreting matter (Cheng and Ruderman 
1989), 

r A « 3 x 10* L^B^ (^j ‘ R\ on cm. (2) 

The electrostatic gap hovers on the accretion disc from the in¬ 
nermost disc radius, ~ r () up to a distance of about r A . Without 
loss of generality we can assume that the stellar magnetic dipo¬ 
lar moment (/i = B.Rl/2) that induces /2-field lines across the 
disc is aligned with the rotation axis of the system (Cheng and 
Ruderman 1991). The //-field on the disc, however, strongly 
depends on the screening factor generated by currents induced 
in the disc surface. In what follows, we take a mean //-field of 
6000 G (Ghosh and Lamb 1979b). 

Protons entering into the gap from the stellar co-rotating re¬ 
gion are accelerated up to energies e V max and directed to the 
accretion disc by the field lines. The maximum current that can 
flow through the gap can be determined from the requirement 
that the azimuthal magnetic field induced by the current does 
not exceed that of the initial magnetic field (Cheng and Ruder¬ 


man 1989) 

7max - 1 .5 x 10 24 /r 2 'Rif L^B^ 1 esu s " 1 . (3) 

The mean number of protons impacting the disc is huge, 
N p = J„rdx/e ~ 10 33 s" 1 , and the total power deposited by the 
proton beam in the disc is P max = y max V max ~ 10 36 erg s -1 . The 
collision of the relativistic proton beam with the disc initiates 
hadronic and electromagnetic showers, in which high energy 
vs are produced from decay of charged pions, whereas 7 -rays 
are produced from the decay of their neutral partners. 

3. PHENOMENOLOGY 

The system just described is an ideal source of high energy 
vs: It has a very dense material for pp—interactions, and, at 
the same time, the region of acceleration is separated from that 
of high density, where interactions occur. The high energy v- 
production in accretion discs of neutron stars is, however, sub¬ 
ject to very stringent conditions on the disc grammage. The 
latter should be large enough as to allow protons to interact, 
and small enough as to allow pions to decay or to escape, so 
as to avoid losing energy. But the disc grammage is actually a 
periodically varying function in these systems, which follows 
the orbital dynamics. As noted, this fact was used to explain 
why X-ray maximum is coincident with a non-detection in the 
7 -ray band for the system A0535+26 (Romero et al. 2001). An 
examination of the cross sections involved is in order. 

In the energy range of interest, the pir cross section, a pn ~ 
25 mb (Carroll et al. 1979), is not too far away from the cross 
section for 7 -ray absorption in the Coulomb field produced by 
disc ions, ~ 10 mb (Cox et al. 2001, Bethe and Heitler 
1934). Then, photons with energies > 100 MeV undergo in¬ 
teractions within typical accretion discs if the mean density of 
hydrogen is n ;> 1 x 10 20 cm -3 (the survival probability for 7 ’s 
to a distance comparable to the thickness of the disc due to 
this process is less than 0.1%). Moreover, assuming a source 
luminosity L37 ~ 1 in the form of ss 1 keV photons and an 
average disc radius of r ss 2 x 10 8 cm, the energy density in¬ 
side the disc is roughly L^/ilnrc) ~ 10 9 erg/cm 3 , imply- 
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ing a number density of 1 keV photons ~ 10 18 cm -3 . There¬ 
fore, the secondary electrons and 7 -rays most effectively inter¬ 
act with the radiation field of the disc. These secondary elec¬ 
trons originating from decays of neutral and charged pi-mesons 
will trigger an electromagnetic cascade in this field resulting 
to a standard cascade spectrum with a cutoff energy around 
m 2 c A /e « 250 (e/lkeV) _ 1 MeV, where m e is the mass of the 
electron and e is the average energy of the background thermal 
photons. 

According to canonical accretion disc models, the average 
density scales as n oc L ^ 35 (M*/M 0 ) -1 / 28 r -21 / 20 (Shapiro and 
Teukolsky 1983). We normalize the particle density of the ac¬ 
cretion disc at 7 in order to fit 7 -ray observations of A0535+26 
(strong photon absorption in the periastron of the system, where 
L37 ~ 1 , Finger et al. 1996). Because of the similarity in the 
cross sections discussed above, whenever 7 -rays are absorbed, 
there will be a much reduced signal in u s. Otherwise stated, the 
hadronic shower will also cool down to sub-100 GeV energies. 
But as the disc density is a function of the position in the orbit, 
before and after periastron passage, there are orbital configura¬ 
tions prone to the emission of both, high energy photons and 
neutrinos. 

To simplify the discussion and get a numerical example of 
the processes involved, we set the half thickness of the disc to 
a constant (h ~ 3.5 x 10 6 cm) and consider an orbital configu¬ 
ration for which L37 ~ 0 . 2 , before or after periastron passage. 
The cross section for pp interactions is ~ 50 mb (Battiston 
et al. 1992). Then, the probability for a pp interaction to take 
place within the disc is 1 —e~ 2h l x f 7> 99%, where X p is the mean 
free path of the relativistic protons. It is worthwhile to point out 
that even though the photon and particle densities are similar, 
the /?7 cross section is two order of magnitude smaller than the 
pp cross section. 

On the other hand, the pion mean free path in this configu¬ 
ration is ~ 10 6 cm, which implies a charged pion survival 
probability —to its decay distance— of P « g-Tcr^/K , w h ere 
T n is the 77 ± lifetime, and F is its Lorentz factor. Then, the prob¬ 
ability for producing u s of energy: 200 GeV and 400 GeV is 
—on average— 10% and 1%, respectively. We expect a rapidly 
fading signal of neutrinos with energies greater than 1 TeV. 

There appears to be enough room for the disc density (or 
grammage) to acquire values such that all protons will partake 
in hadronic interactions, whereas some of the pions produced 
will decay to high energy neutrinos at a level significant enough 
as to show up in neutrino telescopes. To analyze this in detail, 
numerical simulations of a proton-induced cascade unfolding 
within an accretion disc were performed. 

4. NUMERICAL SIMULATIONS 

The accretion disc itself is simulated as a cylindrical volume 
of thickness 2 /z and is filled with elemental hydrogen at a den¬ 
sity «=4x 10 19 cm -3 . A 400 TeV value is retained as the av¬ 
erage energy of each proton in the beam even for L 37 ss 0 . 2 , 
because the uncertainties in all other parameters involved di¬ 
lute the deviation from the fiducial result. The generation and 
tracking of secondary particles in the cascade development was 
performed using GEANT4: a simulation toolkit (designed for 
operation up to center-of-mass energies \fs ~ 100 GeV) that 
provides general-purpose tools for the simulation of the pas¬ 
sage of particles through matter (Agostinelli et al. 2002). We 
process the initial hadronic collisions (with sfs >100 GeV) 
using the event generator DPMJET-II (Ranft 1995). This pro¬ 


gram, based on the Gribov-Regge theory, describes soft parti¬ 
cle interactions by the exchange of one or multiple Pomerons 
(the inelastic reactions are simulated by cutting Pomerons into 
color strings which subsequently fragment into color neutral 
hadrons). In the energy range of interest, \fs ~ 10 2 - 10 3 GeV, 
the average energy fraction of the highest energy baryon in the 
simulated pN collisions is roughly 30%. Therefore, on average, 
after three collisions, the energy of the leading particle will be 
degraded down to ~10 TeV. The secondary meson spectra from 
these three collisions are given in Fig. 2. Additionally, for these 
center-of-mass energies, the soft baryon channel comprises 2 
nucleons, each carrying (on average) 1% of the energy of the 
incoming proton. These secondary distributions together with 
the leading particle —degraded in energy— were injected into 
the cylindrical volume representing the disc. The position of the 
first interaction was selected randomly following a Poisson dis¬ 
tribution with a mean equal to the proton mean free path. The 
second and third interaction points were selected with the same 
procedure, taking the preceding collision as reference, and go¬ 
ing further down within the disc. We then use the GEANT4 to 
track all particles. In this second step, all hadronic collisions are 
processed with GEANT4 implementation of GEISHA (Fesefeldt 
et al. 1985), a program tuned to analyze experimental results 
(of a variety of projectiles and targets) in the few-GeV energy 
range. A magnetic field of 6000 G threads the synthetic disc 
volume, whose direction is pointing upwards in Fig. 1. Vari¬ 
ations in the /i-ficltl across the disc do not affect the results 
presented below, because most of the particles in the shower 
are produced in “minimum bias” events boosted in the forward 
direction (see Fig. 3). 

One can see in Fig. 4 the average number of us produced 
from pion decay per unit energy (bin size = A E v = 50 GeV). 
The tail of the energy distribution is very well fitted by a power- 
law, Ao£“ 7 , with values of Ao and 7 given in the figure. There¬ 
fore, a monoenergetic beam of high energy protons impacting 
onto an accretion disc produces outgoing us with a power law 
spectrum. 

5. THE NEUTRINO SIGNAL 

Having estimated the //-spectrum, we now analyze the 
signal-to-noise (S/N) ratio for a km-scale detector like ICE¬ 
CUBE (Karle et al. 2002). This detector will consist of 4800 
photomultipliers, arranged on 80 strings placed at depths be¬ 
tween 1400 and 2400 m under the South Pole ice. The strings 
will be located in a regular space grid covering a surface area of 
1 km 2 . Each string will have 60 optical modules (OM) spaced 
17 m apart. The number of OMs which have seen at least one 
photon (generated by Cerenkov radiation produced by the muon 
which resulted from the interaction of the incoming u in the 
Earth’s crust) is called the channel multiplicity, N C h. The mul¬ 
tiplicity threshold is set to N c t, = 5, which corresponds to an en¬ 
ergy threshold of 200 GeV (Alvarez-Muniz and Halzen 2002). 
The angular resolution of ICECUBE will be 0.7°, which im¬ 
plies a search window of ~1° radius (Karle et al. 2002). 

The event rate of the atmospheric //-background that will be 
detected in the search bin is given by 

dN f <7$ B . „ 

— 7 - =A e ff / dE I/ ——P I/ ^ tp ,{E v )ATli<y K \o, (4) 
dt B J dE v 

where A e ff is the effective area of the detector, AFli° xl o « 
3 x 10” 4 sr, and d$n/dE v £ 0.2(£„/GeVr 3 ' 21 GeW 1 cm' 2 
s -1 sr~' is the u p + u p atmospheric r/-flux (Volkova 1980, Lipari 
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Fig. 2. — Characteristics of hadronic interactions implemented in the Monte Carlo event generator DPMJET-II. In the first panel we show the average multiplicities 
from pN collisions as a function of the center-of-mass energy (solid lines) together with the energy fraction (in %) going into secondary particles (dashed lines), 
from bottom to top, , tt^ , and 7 . The circles, squares, and triangles give the average numbers of secondary particles produced in non-diffractive pp collisions 
measured by the UA5 Collaboration (Ansorge et al. 1989). The other three panels show the spectra of secondary particles produce in pN collisions at different 
center-of-mass energies. The spectrum of K® L +K® S is roughly equal to that of . 
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Fig. 3.— Neutrino lateral distribution from a proton shower in a typical accretion disc setting. The xy-plane is parallel to the disc. Because the secondary pions 
are produced with very low transverse momentum, most neutrinos are produced near the shower core. 
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Fig. 4.— Neutrino energy distribution generated by pion decay induced by the collision of 400 TeV protons onto an accretion disc of a typical X-ray binary. The 
error bars indicate the RMS fluctuations for each of the mean values, represented by the height of each box in the histogram. The latter were obtained averaging over 
100 showers. The solid line is a fit to this spectrum, whose parameters are shown in the insert. The inset figure is the high energy tail of the distribution, together 
with the same fit. 


1993). Here, P V ^(E V ) « 3.3 x \Q~ n {E v /G€Wf - 2 denotes the 
probability that a v of energy E v ~ 1 —10 3 GeV, on a trajectory 
through the detector, produces a muon (Gaisser et al. 1995). 
On the other hand, the ^-signal is 


= A eff / dE v —^P v ^(E v ) , (5) 

s J dE v 

where 

d* _ MOV 1 1 d<S >Q 

dE v ) 47 rd 2 dE v ( 

is the incoming zx-flux, emitted by a source at a distance d 
from Earth with beaming factor AO/4 tt. Phenomenological 
estimates yield AG/4 tt ~ 0.1 (Cheng and Ruderman 1989). 
The i'-emission spectrum is just d$o/dE u = N p dN/dE v , where 
dN/dE v can be read from the first panel of Fig. 4, and N p « 
3 x 10 33 is the total number of protons impacting onto the disc. 
It is noteworthy that the flux of neutrinos, which is dominantly 
v p + v p at production, is expected to be completely mixed in fla¬ 
vor upon arrival at Earth. Specifically, there is now strong evi¬ 
dence for maximal mixing among all neutrino species (Fukuda 
et al. 1998), which implies that the (/-flux would be completely 
mixed after a propagation distance of (Bilenky et al. 1999) 


dN 

dt 


L 


osc 


2.5 


E v eV 2 
GeV Am 2 


km. 


(7) 


where Am 2 <; 10 6 eV 2 is the zx-mass splitting. Putting all this 
together we obtain 


d<K„ _ 1 d$ 
dE u 3 dE v 

1 We thank Rodin Porrata (LLNL) for bringing this to our attention 


( 8 ) 


A0535+26 is very close to Earth, distant only ~ 2.6 kpc (Gio- 
vannelli and Sabau Graziati 1992). In addition, as only vs go¬ 
ing through the Earth can be identified, the northern location 
of this source makes it accessible from the south pole site of 
ICECUBE. It can be seen in Figure 5 the expected //-flux from 
A0535+26. The atmospheric z/-background for a 1° x 1° bin 
is indicated by the diagonal lines. Also indicated in the figure 
are the fluxes of several others point sources. It is clear that 
for 300 GeV < E v < 1 TeV, vs coming from A0535+26 can 
dominate both over the background and over all other known 
sources. The S/N ratio in the range E v £ [300GeV, 1 TeV], 
using Eqs. (4) and (5) and considering an observing time of 
~ 50 days, compatible with the period in which the accretion 
disc is forming or disappearing, is 


N\ s 

N sM 


1.7. 


(9) 


It is possible, then, that A0535+26 could be detected as a vari¬ 
able neutrino source even within the time span of one orbital 
period. The possibility of detection can be further improved by 
adding up several orbital periods, of which there are three per 
year. It is therefore possible that sources like A0535+26 could 
have been already detected in the AMANDA data. In order to 
search for this possibility two all-sky maps could be used. 1 
One, adding up all observations in periods where the source is 
expected to be ON, before and after periastron passage. An¬ 
other, adding up all observations with the source in OFF mode, 
at the periastron. The difference between these two maps, 
weighted by the respective integration times, should leave only 
one source on the sky, disregarding that this source in itself is 
well below the background for the experiment. Of course, the 
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Fig. 5.— The expected + F /i -flux coming from A0535+26 is represented by the decreasing dash-dotted curve with sharp cutoff above 1 TeV. The flux above 
this energy is subject to large uncertainties because of shower-to-shower fluctuations, and it is not being taken into account to compute the S/N ratio, which is based 
only in the interval (300 GeV - 1 TeV). The shape of the cutoff given in this figure corresponds to that estimated using the phenomenological arguments about pion 
survival probability discussed in the main text. The atmospheric + -background within a 1° x 1° square of the source is indicated by the two diagonal solid 
lines, which corresponds to horizontal and vertical (upper curve) impacts (Volkova 1980, Lipari 1993). Also shown is the current upper limit on flux for a 

source with differential energy spectrum oc E~ 2 and declination angle larger than 40° as reported by the AMANDA Collaboration (Ahrens et al. 2002), as well as the 
expected sensitivity of ICECUBE after its first three years of operation (Hill et al. 2001). Assuming the flux is fairly uniform with decreasing angular bin, we scaled 
down to AQ 2 0 X 2 ° (angular resolution of AMANDA) the current upper limit on diffuse -l-F^-fluxes (90% CL) as reported by the Frejus Collaboration (Rhode et 
al. 1996). We stress that the Frejus detector was located in an underground laboratory in the Alps and consequently it had not a complete field of view of A0535+26. 
The decreasing dash-dotted line indicates the - spectrum of SN IC443 (Hill et al. 2001). This curve relies on the assumption that (on average) a pair of muon 

neutrinos with energy E n /4 are produced in the pion-to-muon-to-electron decay chain of proton showers (in very low density media) for every photon with energy 
E n /2 (Gaisser et al. 1998). The bumpy curve indicates one of the hypotesyses (Model-II in Bednarek and Protheroe 1997) for -t-F^-production in the Crab 
Nebula. 
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same technique can work for ICECUBE. A few years integra¬ 
tion time would add up tens of orbital revolutions, and detection 
of A0535+26 should be unambiguous. Indeed, A0535+26 is, in 
projection, less than 5° away from the Crab Nebula, which was 
suggested as a potential detection by AMANDA (Barwick et 
al. 2002). This detection is based on a ~ 6 ° x 6 ° binning of the 
nothern hemisphere. It is not implausible, then, that A0535+26 
may significantly contribute to this angular bin, particularly at 
energies below 1 TeV. 

From a technical point of view, ICECUBE must select a soft 
energy cut-off to observe this source, since A0535+26 shines 
mostly below 1 TeV (the channel mutiplicity threshold of the 
experiment must be set as 5 < (V c h < 30 (Karle et al. 2002). 
A next generation observatory, like NEMO (Riccobene et al. 
2002), would improve the angular resolution to ~ 0.3°, which 
would increase the .S'/A-ratio by a factor of three. 

6 . CONCLUDING REMARKS 

In summary, a wide range of parameters able to describe 
the physical situation of accreting neutron stars make them a 
bona-fide population of ^-emitters with power-law spectra. In 
the particular of case of A0535+26, the + f^-flux between 
300 GeV and 1 TeV overwhelm those of all other point sources, 
and is clearly above the atmospheric background in a 1 ° x 1 ° 
angular bin. The neutrino signal will be periodic in nature. 


Alvarez-Muniz J. and Halzen F. 2002, Astrophys. J. 576, L33. 

Agostinelli S. et al. [GEANT4 Collaboration] 2002, SLAC-PUB-9350. 

http://wwwinfo.cern.ch/asd/geant4/geant4.html 
Ahrens J. et al. [The AMANDA Collaboration] 2002, arXiv:astro-ph/0208006. 
Ansorge R.E. et al. [UA5 Collaboration], Nucl. Phys. B328, 36. 

Battiston R. et al. [UA4 Collaboration] 1982, Phys. Lett. B 117, 126. 

Barwick S. W. et al. [AMANDA Collaboration] 2002, arXiv:astro-ph/0211269. 
Bednarek W. and Protheroe R.J. 1997, Phys. Rev. Lett. 79, 2616. 

Berezinsky V.S., Bulanov S.V.. Dogiel V.A., Ginzburg V.L., and Ptuskin V.S. 

1990. Astrophysics Of Cosmic Rays, North-Holland. 

Berezinsky V.S., Castagnoli C. and Galeotti P. 1985, Nuovo Cim. 8C, 185 
[Addendum-ibid. 8C, 602, 1985]. 

Bethe H.A. and Heitler W. 1934, Proc. Roy. Soc. London A146, 83. 

Bilenky S.M.. Giunti C. and Grimus W. 1999, Prog. Part. Nucl. Phys. 43, 1. 
Carroll A.S. et al. 1979, Phys. Lett. B80, 423. 

Cheng K.S. and Ruderman M. 1989, Astrophys. J. 337, L77. 

Cheng K.S. and Ruderman M. 1991, Astrophys. J. 373, 187. 

A. N. Cox (Editor) 2001, Astrophysical Quantities , Fourth Edition, Springer- 
Verlag, New York, p.213-214. 

H. Fesefeldt 1985, PITHA-85-02, CERN-DD-EE-81-1, CERN-DD-EE-80-2. 
Finger M.H., Wilson R.B. and Hannon B.A. 1996, Astrophys. J. 459, 288. 

Y. Fukuda et al. [Super-Kamiokande Collaboration] 1998, Phys. Rev. Lett. 81, 
1562. 

Gaisser T.K. and Stanev T. 1985, Phys. Rev. Lett. 54, 2265. 


The accumulated signal during the apoastron -where there is 
no disc formed- and periastron passages will always be at the 
noise level. On the contrary, the accumulated signal during 
time spans before and after the periastron passage, will be, as 
we have shown, significantly different from noise level. A rea¬ 
sonable integration time for neutrinos would then secure high- 
confidence detections of the appearance and disappearance of 
the accretion disc (> 3yr to reach a 5a effect). Complementary 
information from high-energy photon astronomy in the X- and 
7 -ray domain will be essential to confirm the hadronic origin 
of the radiation in accreting neutron star systems, as well as 
to study the formation and loss of accretion discs in eccentric 
binary systems. 
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